Collective Oscillations of Strongly Correlated One-Dimensional Bosons on a Lattice 
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We study the dipole oscillations of strongly correlated ID bosons, in the hard-core limit, on a 
lattice, by an exact numerical approach. We show that far from the regime where a Mott insulator 
appears in the system, damping is always present and increases for larger initial displacements of the 
trap, causing dramatic changes in the momentum distribution, nu- When a Mott insulator sets in 
the middle of the trap, the center of mass barely moves after an initial displacement, and nu remains 
very similar to the one in the ground state. We also study changes introduced by the damping in 
the natural orbital occupations, and the revival of the center of mass oscillations after long times. 
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Experimental realization of one-dimensional (ID) 
atomic gases has opened the possibility of studying 
equilibrium and non-equilibrium properties of systems in 
which quantum fluctuations play a fundamental role. For 
example, the ID superfluid-Mott-insulator transition 
was observed experimentally by Stoferle et al. In 
addition, ID geometries have allowed the observation of 
the strongly correlated hard-core boson limit with and 
without ja] a lattice along the ID tubes. 

Recent experiments studying transport properties of 
trapped ID Bose gases 00 have reported a large damp- 
ing of the center-of-mass (cm.) oscillation when an axial 
lattice is present in the system. Additionally, an over- 
damping was observed in which the cm. remains dis- 
placed from the middle of the trap 0, • In Ref- 0] this 
effect was related to the presence of an incompressible 
Mott insulating phase, while Fertig et al. reported 
to be far from this regime. For weak interactions, and 
small trap displacements, the damping has been argued 
(within mean field approximations) to relate to the exis- 
tence of non-condensate fractions, which cause a "dissi- 
pative" behavior of the cm. oscillations B, El. For large 
trap displacements, dynamical instabilities |lft IrH ] can 
also produce damping. 

For strongly correlated ID bosons on a lattice, the 
mean-field description assuming Bose-Einstein Conden- 
sation (BEC) breaks down. For hard-core bosons (HCB), 
the largest eigenvalues of the one-particle density matrix 
(OPDM) [also called natural orbitals (NO)] scale pro- 
portional to \/Wb [Nb is the number of bosons]. There 
is no BEC. However, there is a finite supcrfluid fraction 
p s , which for the periodic case (no trap) can be easily 

in sm( tv N b /N) 

calculated pi = ; , where N is the total 

nN b /N 

number of lattice sites (see e.g. Ref. |l2j). In two di- 
mensions (2D) HCB exhibit BEC at zero temperature. 
However, the 2D superfluid fraction is smaller, at any 
density, than in ID, though larger than the 2D conden- 
sate fraction 01 ■ These results exemplify the non-trivial 
differences between BEC and superfluidity. It is unclear 
how much the lack of a true BEC is central to under- 
standing transport properties in ID. 



The presence of the trap in the experiments with ultra- 
cold quantum gases generates an inhomogeneous density 
profile 0, Q| , and modifies the single particle spectrum 
of the periodic lattice ^{|. This could, in general, change 
the behavior of the superfluid and condensate fractions. 
However, for the ID HCB, the power-law decay of the 
OPDM (p x ^ 1/ ,/x) , and the occupation of the lowest 
NO (A ~ \fNb~), have been shown to be the same as in 
the periodic system [lij . 

In this work we study the dynamics of the cm. oscil- 
lations of hard-core bosons (HCB) confined in ID har- 
monic traps with an underlying lattice, by an exact nu- 
merical method. Our purpose is to examine the damping 
of such oscillations, and its consequences in the momen- 
tum distribution function (rik) and natural orbital (NO) 
occupation. We show that far from the Mott-insulating 
(MI) regime in the middle of the trap, the damping of 
the cm. oscillations grows with initial displacement and 
in the large damping regime, causes dramatic changes in 
rife. On the contrary, when the MI appears in the middle 
of the trap, even for small initial displacements the cm. 
barely moves from its initial position, and remains 
very similar to the one in the ground state. 

The HCB Hamiltonian can be written as 



H = —t 



bib. 



h.c. 



(1) 



with the addition of the on-site constraints b\ 2 = b\ = 0, 

= 1. These constraints on the creation (b\) and 

annihilation (bi) operators avoid double or higher occu- 
pancy. In Eq. QJ, the hopping parameter is denoted by t, 
and the last term describes a harmonic confining poten- 
tial with curvature V2. n» = b\bi is the particle number 
operator. The non-equilibrium dynamics of the system, 
when at time r = the trap is displaced a distance xo, 
is obtained by means of the exact approach introduced 
in Ref. 0|. It is based on the Jordan- Wigner trans- 
formation, and allows obtaining the HCB one-particle 
Green's function from which all the quantities analyzed 
in this work are calculated. In addition to being exact, 
this method allows studying relatively large system sizes. 
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Previous studies have shown that the key parameter 
that controls the thermodynamic behavior of this model 
is the characteristic density p = N h ^fV^Jt Q. As p — > 0, 
one recovers the continuum limit . On the other hand, as 
p increases beyond a critical value (p c ~ 2.6 — 2.7) a MI 
region builds up in the middle of the trap, where the 
density is pinned to unity with no quantum fluctuations. 

In Fig. QJa)-(c) we show the cm. oscillations of a 
trapped HCB system for p = 1 < p c , far from the regime 
with a MI in the middle of the trap. At time t = the 
maximum density in the trap is n = 0.48, and the cloud 
radius is Rq ~ 145a, with a the lattice constant. Fig.QJa) 
shows that even for a very small displacement, 7% of the 
cloud radius, a damping of the oscillations occurs. For 
this displacement the maximum momentum of the cm. 



of the system is ~ 7r/30a, very far from the momentum at 
which the classical modulation instability occurs, which 
for HCB is ir/2a (see e.g. Ref. [HI). Since HCB can be 
mapped into noninteracting fcrmions, one can immedi- 
ately see that damping occurs for any initial density and 
trap displacement. In the fermion language, the damp- 
ing is caused by a dephasing of the particles, which arises 
due to their dispersion in the lattice-harmonic trap sys- 
tem 0, 0] . As the initial displacement of the trap is in- 
creased [Fig.^b)], the damping also increases. For large 
displacements, Bragg scattering starts to occur produc- 
ing a shift of the cm. oscillations [Fig. dc)]. A large 
damping of the cm. oscillations can also be observed. In 
Fig.^c) the maximum momentum of the cm. is ~ 7r/4a, 
still smaller than it /2a. 
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FIG. 1: (color online). Evolution of the cm. position (a)-(c), maximum value of nt (d)-(f), and the occupations of the three 
lowest NO (g)-(i), vs time (r, in units of h/t). The initial trap displacements are xo = 10a [xo/Ro ~ 0.07] (a),(d),(g), xo = 29a 
[xo/Ro ~ 0.2] (b),(e),(h), and x = 80a [x /Ro ~ 0.55] (c),(f),(i). These systems have N b = 101, and V 2 a 2 = 10" 4 t. 



We next analyze the consequences of the damp- 
ing of the cm. motion in another physical ob- 
servable, the momentum distribution function = 
(a/OEji^^PM, with C = (V 2 /t)- 1/2 0, and 
Pji = (b]k) the OPDM. In Fig. Gld)-(f) we show how 
the maximum value of (n™) evolves as a function of 
time. In Fig.^d) one can see that even for small damp- 
ing, the oscillations of n™ are accompanied by an overall 
decrease of its mean value, which reaches almost half of 
the original height for r = 4000fi/£. The changes in n™ 
become more dramatic with increasing xq. In the large 
damping regime [Fig. ^f )] , one can see that n™ reduces 
almost to its minimum value in the first oscillation of the 



cm. The reduction of is accompanied by a large in- 
crement in the full width at half maximum w, as shown 
in Fig. |21 These results are in contrast with the ones for 
the equivalent noninteracting fermions [l5j . where after 
the damping of the cm. oscillations, n™ and w are al- 
most the same as the ones at r = 0. (Contrary to the 
density, the evolution of rife for fermions and bosons is 
very different.) 

At this point it is relevant to know the behavior of the 
"condensate" occupation during the oscillations of the 
cm. For these systems it is important not to confuse 
with the "condensate" occupation. The condensate can 
be defined as the largest eigenvalue of the OPDM (the 
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lowest NO) [HII3, 

and since for HCB its occupation Ao 
scales proportionally to y/Nb we prefer to call it a 

quasicondensate (Ao — ► oo for Nt — ► oo, but Xo/Nf, — > 0). 
One particular example in which the differences between 
rtl" and Ao are extreme was reported by one of us in Ref. 

. During the expansion of the HCB gas, after turning 
the harmonic trap off, decreases while Aq increases. 
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FIG. 2: (color online), nt of the ground state (GS) compared 
with the ones at r = 4000h/t for two different initial displace- 
ments of the trap (Nb = 101, and Via 2 = 10~ 4 t). As shown in 
Figs.|H b ).( c )> at r = 4000h/t the motion of the CM is totally 
damped. The inset shows the OPDM for the same systems, 
the straight line is pij 



,-l/a. 



In Fig.njg)-(i) we show the time evolution of the three 
lowest NO occupations (the highest occupied ones). It 
can be seen that indeed the behavior of Ao is different 
from the one of n™. At short times, and for low damp- 
ing rates [Fig. mg)], the lowest NO occupation does not 
change during the oscillations of the system, even though 
the lowest NO itself does oscillate in a coherent man- 
ner, as would be the case for an ideal BEC when there 
is no lattice. Only when the cm. energy is transfercd 
into excitation energy, the lowest NO occupations (and 
the shape of their wave-functions) start to change. This 
reflects the onset of a cut-off to the off-diagonal quasi- 
long range order initially present in the OPDM, and it 
is similar to the appearance of a finite temperature in 
the system. Such effect can be seen in the inset in Fig. 
121 where we show the OPDM of the ground state when 
the cm. oscillations are damped. In the latter cases, the 
fast decays of the OPDM reflected in the plots are very 
close to exponentials. As expected, the largest damping 
produces the fastest decay. The consequent reduction of 
the lowest NO occupations can be seen in Fig. IHh),(i). 

So far we have analyzed trapped systems with no MI 
phase. In what follows we study the consequences of 
approaching the MI regime on the damping of the cm. 
oscillations. Fig. |3fa) shows that after damping, when 
the characteristic density of the system is increased, the 
final cm. position is displaced from the center of the 
trap. We kept the ratio between the initial trap displace- 
ment xq and the cloud radius Rq constant, and equal to 
0.2 as in Figs. 2Ib),(e),(h). The insets show the results 
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FIG. 3: (color online), (a) Ratio between the cm. position 
after damping (xf) and the initial trap displacement (xq) vs 
p. The insets show xcm /Ro vs r for the signaled p. (b) Ratio 

=0 ) and n fe = at r = (n° =0 ) 



between 7ifc=o after damping (n k - 



vs p. The inset shows n™ vs r = for the signaled p. 

for xcm/%0 vs t hi traps with p = 2.5, just before the 
MI phase appears, and for p = 4, when a MI domain 
is present in the center of the trap [Fig. E[a)]. These 
results show that (i) the damping rate increases with in- 
creasing p, and Bragg scattering keeps the cm. displaced 
from the middle of the trap, (ii) when there is a MI do- 
main (p = 4), the cm. barely displaces from its original 
position. These results are in agreement with the ex- 
periments reported by Stoferlc et. al Q, and differ from 
the ones in the weakly interacting regime where for small 
initial trap displacements no shift of the cm. position is 
observed after damping 0,0. 
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FIG. 4: (color online). Density and momentum profiles for a 
system with Nb = 101, and p = 4. See also results for p = 4 
in the insets of Fig. [3] 

One interesting feature that appears with the forma- 
tion of the MI can be seen in the inset in Fig. UJb). 
Although n™ at r = is small compared with njj 1 for 
smaller p, its value remains almost the same during the 
evolution of the system. The same occurs with w, the 
width of rifc, as seen in Fig. Ub). The above result can 
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be intuitively understood by the fact that almost no cm. 
energy has been converted into excitation energy, which 
keeps the system in a state with n k , and NO occupations, 
similar to the ones in the ground state. Something simi- 
lar may occur in experiments in the overdamped regime, 
when the cm. almost does not change from its initial 
position 0- However, we should remark that we do ob- 
tain big changes of n™, and w, [Figs. El and El^b)] as 
a consequence of the damping in systems with no MI. 
This is in contrast with the behavior of w one infers from 
the constant cloud widths after time-of-flight (TOF) re- 
ported experimentally in Ref. 0. The reason for this 
difference may be that TOF measurements, due to inter- 
particle interactions during the expansion, may not give 
the original n k in the trap. (An explicit example was 
presented in Ref. If this is the case, other experi- 

mental techniques will be needed to determine the effects 
of the damping on rifc . 
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FIG. 5: (color online). Long time evolution of the cm. posi- 
tion (a), maximum value of nu (b), and the occupation of the 
lowest NO (c), vs time (r, in units of h/t). The initial trap 
displacement is xo = 4a [xo/Ro ~ 0.1], and the system has 
N b = 20, and p = 1.0. 

One final remark is in order concerning the damping 
discussed here. Since the experimental systems with ul- 
tracold quantum gases are almost closed systems, damp- 
ing has a different meaning when compared with systems 
in contact with a reservoir, where the energy is truly dis- 
sipated. In the systems we have analyzed, the energy is 
conserved. In Fig. we show a long time evolution of 
the cm. oscillations, n™, and the lowest NO occupation, 
after an initial displacement of the center of the trap. As 
seen in these figures, a revival of the cm. oscillations 
occurs, with an increase of nj™ and Ao to values similar 
to the ones at r = 0. This suggests that ultracold gases 
experiments could be used to study quantum Poincare re- 
currences |2^ | . The relation between long time dynamics 
in integrable systems, such as HCB analyzed here, and 



nonintegrablc (finite on-site repulsion) models, is an im- 
portant open question (of broad relevance to condensed 
matter physics) that needs further analysis. 

In summary, we have presented a detailed, and exact, 
study of the damping of the cm. oscillation in a sys- 
tem of strongly correlated bosons. We have shown that 
far from the Mott-insulating regime large damping oc- 
curs only for large initial displacements of the trap, and 
it reduces dramatically the maximum value of rik and 
increases its full width. The damping also destroys the 
quasi-long range order of the OPDM, reducing the oc- 
cupation of the lowest NO, analogous to an increase in 
the temperature of the system. When the Mott insula- 
tor is present in the trap, overdamping occurs even for 
small initial displacements. In this case the cm. remains 
very close to its initial position, and no big changes occur 
in rife. Finally, we have shown that since the system is 
closed, after some time a revival of the cm. oscillations 
occurs, and n™ and Ao also regain their initial values. 

Note added. Complementary to this work, a related 
study has been recently presented in Ref. |24| . 
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